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Fifth Quarterly Report 

STUDY ON SPECTRAL/RADIOMETRIC CHARACTERISTICS 
OF THE THEMATIC MAPPER FOR LAND USE APPLICATIONS 


1 . OBJECTIVE 

The objective of this investigation is to quantify the performance 
of the TM as manifested by the quality of its image data, in order to 
suggest improvements in data production and to assess the effects of 
the data quality on its utility for land resources applications. Three 
categories of this analysis are: a) radiometric effects, b) spatial 
effects and c) geometric effects, with emphasis on radiometric effects, 

2. TASKS 

Four tasks have been established to address the above objective. 

The first three are to study radometric performance, spatial performance 
and geometric performance, respectively , while the fourth is to study 
spectral characteristics. In keeping with the identified objective, 
the radiometric performance study is the major task. 

3. STATUS AND TECHNICAL PROGRESS 

During this fifth quarterly reporting period, work was continued on 
our study of the radiometric performance of TM. Additional characteri- 
zation of scan-related level -shift noise in TM reflective-band data was 
achieved through analyses of non-thermal data from night scenes. Two 
other analysis efforts utilizing the nighttime data were initiated. The 
first is a resumption of our earlier analysis of scan-direction-related 
signal droop effects and the second is an examination of high-frequency 
noise within scan lines. Also we responded to a request for inputs by 
making recommendations for the procedure used in TIPS (Thematic Mapper 
Image Processing System) to radiometrically calibrate TM data prior to 
their distribution. Also, technical presentations were made at the 
Pecora VIII Symposium and the December LIDQA Investigators' Workshop. 

A paper was prepared and accepted for publication in the IEEE Trans- 
actions on Geoscience and Remote Sensing. 

3.1 PROBLEMS 


None. 
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3.2 ACCOMPLISHMENTS 

Accomplishments In four technical areas are described below. 

3.2.1 ANALYSIS OF LOW-FREQUENCY LEVEL-SHIFT NOISE 

As reported previously, examination of nighttime reflective-band 
data from the Augusta night scene proved very fruitful in analysis of 
low-frequency level-shift noise. Analysis of the Buffalo night scene 
during this quarter provided further insight into the characteristics 
of that noise. In this scene, the Form #1 noise (most evident in 
Band 1 Detectors 4, 12, 10, and 8) is absent, while the Form #2 noise 
(exemplified by Band 7 Detector 7) is present to the same extent ob- 
served in other scenes. In the seven scenes we have examined, Form #2 
noise is always present, and Form #1 is present in most, but not all, 
scenes (See Table 1). The peak-to-peak amplitude for each detector 
for each form of the noise is essentially constant in all cases where 
that form of the noise exists. The phase relationships between the 
affected detectors also remain constant in all images (i.e., Band 7 
Detector 7 is always in its "high" state when Band 5 Detector 8 is in 
its "low" state). Tables 2, 3, and 4 list the magnitude and phase of 
each form of the noise for Scenes 40037-02243, 40049-16262, and 
40161-02481, respectively. The poorer separation of states shown 
for Frame 40049-16262 is due to the fact that these tables were made 
from scene radiometric data, and this scene was the only daytime scene 
of the three. 

3.2.2 ANALYSIS OF SCAN-DIRECTION* RELATED SIGNAL DROOP EFFECT 

Initial analysis of the nighttime reflective data for signal droop 
effects is promising. Although the magnitude of the droop is quite, 
small at these low signal levels, it is clearly and consistently there. 

Further efforts directed at characterizing the effect and proposing 
a solution are planned. 

3.2.3 ANALYSIS OF HIGH-FREQUENCY NOISE 

The success in quantifying the low-frequency noise through the use 
of night data suggested its potential value in analysis of the high- 
frequency noise reported earlier by other investigators. Analysis of 
fast Fourier transform data is still underway, with preliminary results 
supporting previous findings of components with wavelengths of <3.24 
pixel widths (32.1 KHz) and 13.67 pixel widths (7.61 KHz). 
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3.2.4 RADIOMETRIC CALIBRATION PROCEDURE FOR TIPS PROCESSING 
OF TM DATA 

One step In the calculation of Radiometric Look-up Tables (RLUT's) 
for converting received sensor signals into digital counts on radio- 
metrically corrected computer-compatible digital tapes (CCT‘s), utilizes 
absolute calibration information. As the sensor response changes, the 
number of radiometric units per raw data count changes. However, the 
desire is to have a consistent number of radiometric units per CCT 
digital count in each band. A special meeting was called by Dr. John 
Barker during the December 1983 LIDQA Investigators' Workshop*, W. Mai 11a 
and M. Metzler of ERIM attended the meeting. Dr. Barker briefly reviewed 
‘ for the assembled group the history of maximum radiance (R max ) values 
that (a) had been used ii SCROUNGE processing of TM data during the 
first year, (b) have been used in subsequent processing by TIPS, the 
Thematic Mapper Image Processing System, (c) are proposed for the next 
update of TIPS, and (d) may be used for Landsat-5 data. He asked for 
comments on the suitability of the proposed values (or approach for 
thei" selection) and any other suggestions. 

Because of conflicting user requirements, thire is not a clear-cut 
choice. Figure 1 hel^s illustrate the problem * through lines repre- 
senting use of two different R max values fr v he linear relationships 
currently in use. (Rmin values equal to zer».« are assumed for simplicity 
in the example.) The line with the larger value covers a greater 
range of signal radiance values but, because the number of digital counts 
remains the same, produces a coarser quantization of radiances (more 
radiance units per digital count). A large radiance range is desired 
by those researchers interested in mapping snow, clouds, and other 
bright scene materials and by Landsat systems people who wish to allow 
for possible different responsivity in the Landsat-5 Thematic Mapper. On 
the other hand, a finer quantization of radiances is desired by the 
majority of investigators, who are not concerned with accurate repre- 
sentation of radiances from bright objects, and particularly by those 
investigators with bathymetry and other applications for which low 
signal values are of prime interest. 

Dr. W. Mai i 1 a made a suggestion in the meeting that consideration 
be given to use of non-linear (perhaps piecewise-1 inear) calibration 
curves rather than the simple linear relationships that have been 
employed up to this point in time. Figure 2 illustrates a two-piece 
curve that has the following features: (a) radiances up to the rela- 

tively high R ffla x2 value can be represented and (b) the finer quantiza- 
tion associated with R max -| can apply to the majority of data values. 

The sacrifices are a large radiance quantization unit for the high 
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range of signal values (larger even than that associated with Rniax2 
and the simple linear relationship of Figure 1) and a loss of the simple 
linear relationship. To specify the curve, only the two end points and 
the break point need be Identified, e.g., (S n ,L n ), {$, ,L-.), and (So,L 0 ) 
on Figure 2(b). u 0 22 


In discussing the meeting with Mr. Fred Tanis and others at ERIM 
who are involved in bathymetry studies, the desirability for a three- 
piece curve having a modified low-radiance response characteristic was 
identified. The illustration in Figure 3(a) has the shallowest slope 
(most responsive, i.e., most counts per radiance unit) for low radiance 
values, the moderate slope (comparable to the original simple linear 
relationship) for mid-range values, and the steepest slope (least 
responsive, i.e., fewest counts per radiance unit) for high radiance 
values. This type of response characteristic can be specified by 
giving the two end points and the two intermediate break points, as 
shown in Figure 3(b). 


If two- or three-piece characteristics are adopted for generating 
the radiometric lookup tables, the appropriate end points and break 
points should be identified in the headers of CCT's produced by TIPS. 
One other piece of Information should also be included if any detectors 
are more responsive than the selected response characteristics for the 
band. This added information would be the signal count level (after 
correction) which corresponds to saturation of the raw detector output. 


3.3 SIGNIFICANT RESULTS 


See Sections 3.2 and 3.5. 

3.4 PUBLICATIONS AND PRESENTATIONS 

A paper was presented at the Pecora VIII Symposium held in Sioux 
Falls, SD, on October 4-7, 1983, in the session on Landsat 4 results. 
Entitled "Radiometric Analyses of Landsat-4 Digital Image Data", by 
William Malila, Daniel Rice and Michael Metzler, the paper was presented 
by W. Malila. A written paper was submitted for publication in the 
symposium proceedings. 

William Malila and Michael Metzler attended the LIDQA Investigators' 
Workshop at NASA/GSFC on 6-7 Dec 83. W, Malila presented a summary of 
ERIM's now-completed MSS investigation for LIDQA (Contract NAS5-27254), 
while M. Metzler presented an update of our TM efforts under this 
contract. 
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A paper describing results of these investigations was written 
and accepted for publication in the May 1984 issue of the IEEE Trans- 
actions on Geoscience and Remote Sensing. Entitled "Characterization 
of Landsat-4 MSS and TM Digital Image Data", it was co-authored by 
William A. Malila, Michael D. Metzler, Daniel P. Rice, and Eric P. Crist. 

3.5 RECOMMENDATIONS 

As discussed in Section 3.2.4, it is recommended that NASA imple- 
ment a two- or three-piece (the latter being preferable) piecewlse- 
linear calibration characteristic for radiometrically calibrating TM 
data In TM ground segment (TIPS). This would avoid unnecessary satura- 
tion of output values and enhance low-radiance responses, while pro- 
viding for flexibility for accommodating Landsat-5 TM data. 

3.6 FUNDS EXPENDED 

A total of approximately $20,600 was expended during the three 
months September through November 1983. The cumulative spending 
through November represents approximately 54% of the total contract 
award and 77% of the funds allocated. Expenditures during the period 
1-20 December 1983 are not included in these percentages. 
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FIGURE 1 . ILLUSTRATIONS OF RADIOMETRIC CALIBRATION USING SIMPLE LINEAR 
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FIGURE 3. ILLUSTRATIONS OF THREE-PIECE PIECEWISE-LINEAR RESPONSE CHARACTERISTIC 
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Table 1. 

Presence 

of Low- 

-Frequency Level 

Shifts 

in TM Data 

Location 

Date 


Scene ID 

Form 1 

Form 2 

Arkansas 

22 Aug 


40037-16034 

yes 

yes 


Buffalo 

22 Aug 

(PM) 

40037-02243 

no 

yes 


Iowa 

3 Sep 


40049-16262 

yes 

yes 


N. Carolina 

24 Sep 


40070-15084 

yes 

yes 


Cape Cod 

8 Dec 


40145-14492 

yes 

yes 


Georgia 

24 Dec 

(PM) 

40161-02481 

yes 

yes j 

periodic 

Grand Bahamas 

14 Jan 


40182-15125 

no 

yes ( 

periodic 


Table 2. Magnitude of Level Shifts for Might Scene 40037-02243 
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Negative amplitudes indicate level shifts with phase shifts of 180 deg. relative to B1D4 (Fora 1) or Bf£»7 (Fora 2). 

ERIM 6 Dec 83 


Table 3. Magnitude of Level Shifts for Scene 40049-16262 


ORIGINAL PAGE 
OF POOR QUALITY 


p)ww(s«fin»-’-u:CiW^into«3 


* * * » I »»>«•»*»•« ♦ 

r* ’■,» r* t- r* r* t*> Cl 1 * IN T* ^ f“ +* 


nMDNcoMiO'*(oin(<jo)[g<o(DO) Nto«n<oo)«j-ocJiON(sr)^0'r-N '-'•nnrino'riflVNffl^w'-M 

« i i > i i * • t » i ■ * i * > i t > i i i > t * i i i » * * » 

r» r- *■*" P) *“*“»■* P) r* 


InPiVMipininp)o>pJip*-nr*»tncN 
0"0000 O'- 00600000 

ddoddo d d do ddoddo 

1 1 1 1 1 1 1 1 1 1 1 tit 


nuffflVMINWOllIJNMDMMM 

lDOhtfnMNP)^(D(OV’-01W** 

o*-oooo”wor'opQO*-*- 

6 6 ddoddo od ddddod 

1 1 < 1 1 1 1111 


iDr>m<vc<(D(o<Nr>rno>r>r-r'*iDa> 

OOOflOBJItNIIIHONONN 

t 1 1 t 1 1 t r 1 r 1 1 1 < t 1 

OOOOOOOOOOOOOOOO 

I I I I I I I I 


•j'int , )oi«Mtnoow»-ooMN'‘ 

u)r*uj^iNtn!>’”»cvr'OmcDvo 

Ct(N»“0«N''-0i-r5CMr'»WC)Pn'> 

o o d o d d o o o o 6 6 6 6 6 6 


cotooco't in Min’- in ooaronin*- 

C’WtflNOON’-WNWlOUJOOVO’ 

0*'’-*-^000*- 00 0*-OCv|© 

ddoddo o’ ddddodd do 

1 1 1 1 1 1 1 1 1 1 


coM’-toaoxj-Ncnxrromtoeovioifl 

Nwr-*r>OT''r-(flcnoini-t>cocs|co 

ooooo«-ooo*-ono’-o*- 

o o’ o 6 6 6 6 6 6 6 6 6 6 6 6 6 

11 1 1 1 1 1 


** cs co <»• in co r»- uo o> © *- w p> v in to <- w o v in 0 r< w o> o *- c\i p> ** in in *- n cj v in to r- eo o> o r- cv co «r m io - 

T" W S!? T= 5 = *- r« ¥• V"* W V“ 

■<t «t- v >» «r ^ *» * * 't **»■ v ^ *t in in in in in in in ib in in in in in in in m r* r- r- f~ r- r t> r» r- r- r» r- r* r- r> 0 


(niDinNr'<Ncov»“wr3'N*-«j-nin >-nin(N'*>ifO)nf'ifloi)NO<fO*- (ooorotnio^tst^inMrov — ’-w 


CN " 


1- »- V- CN IN t'l «- T“ 


oifVciflONffluioinnioiDiD- tocniDcowisv-eo'Viov-Mcnncoo cotointNOr'OinNeor-Nv-rs.^t^ 
r» cv i fl - r*' r " ** »'*-ct ,, '' r ”tN'^’“ i- ^ cs ^ cn cn 


incnioco’-Of'CiNintatnPoiflwtn 

floiowiniotrcococowr.row^-coo 

'i-T-^Or-OOOOOOOO^Or- 


oin'i-aiOo)io^in«t , <tPotocor»>t 

omvcocoin'j'incnr-vpoininiDLn 

COONO^Ot-OOO'-O’-O’-'- 


mcoc>tinr«iococoo»“iDO<ONtnc^ 

r-tocs«3-ino)^rm»-tj)r-co'-ocoo 

uinnoNOwO’-O’-NMNom 


OOOOOOOOOOOOOOOO 

III I I I I I I I I I I I 


OOOOOOOOOOOOOOOO 

I I I I I I I I 


OOOOOOOOOOOOOOOO 

I I I III I I I 


inCD«ini-ncpc>jr-*t’d , T*ineu<* , j!'> 

r.i~ot>»-r'r'r)rjinw'!rN(n(j)N 

WO’-01'-OOUJT-(I)rMi)--r-O 


M(DO)T-otnco«« , eor-tD'i-coinDo 

NifliaNrsO’-ninwowncinv 

fl>-0fl0’-0l"0000'-rr r 


000 " 0000000"0000 


OOOOOOOOOOOOOOOO 

I 1 I 


eo<f«t-NoiO)f'Oin'tcM<oiNij-in 

otai^rowr-oor-d’-cnO’-ocor- 

inNNOnN’-oonoNnnn'i 

dddddddddddddddo 


»-cMOxtiniflNcoo)o»-NM«-iniD »-«o'tinioiMOoiO^-'iW'twu) ftso'U'intDr'cooiO’-cvo'j'inio 


^. T . r . r --r. r . r . r . r . r . T .„ r . T . r „ r . (SNNfNNISNNNNCSPINMNW rjWOMOOPJOOnjnOOOPJP) 


10 


Table 4. Magnitude of Level Shifts for Might Scene 40161-02481 
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** Negative amplitudes indicate level shifts with phase shifts off 180 deg. relative to B104 (Form 1> or B7D7 (Frna 2). 
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